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1. INTRODUCTION 


The need for cost-effective and reliable space transportation systems has been accentu- 
ated with increasing world-wide competition to exploit space applications and the urgency 
of military payloads. The use of aerodynamic rather than propulsive forces to perform 
various types of orbit transfer can provide significant advantages in increased payload, re- 
duced fuel weight and reduced launch activities. The economic superiority of aeroassisted 
orbital transfer over all-propulsive maneuvering has been demonstrated through various 
studies [1]-[13]. 

In the typical aeroassisted coplanar orbital transfer shown in Figure 1, the vehicle is 
initially in high earth orbit at radius, rj, when a retro impulse, say A Vi, brings the vehicle 
into an elliptical orbit with perigee, say r 9 . Whereas ideally, r p would be the top of the 
atmosphere, for practical reasons it would be selected somewhere inside the atmosphere in 
a realistic aeroassisted orbital transfer maneuver. 

The portion of the orbital transfer of particular interest in this study is the atmospheric 
flight portion which starts as the vehicle enters the atmosphere at a radial distance from 
the center of the earth, say r a . By appropriate modulation of the lift and drag forces, 
the vehicle reduces its speed to a level corresponding to its final lower orbit radius, r 2 . 
The atmospheric trajectory also determines the heating rate which the vehicle skin will be 
subjected to. Thus, it is necessary to select a trajectory which does not subject the vehicle 
to temperature levels higher than can reasonably be accommodated by the vehicle skin. 
At atmospheric exit, the vehicle enters an elliptical orbit with apogee at rj. A circularizing 
impulse, AV*, at the apogee puts the vehicle in the desired final orbit. 

Aeroassisted coplanar orbit transfer has been studied in [4] where it is determined that 
a zero flight path angle at atmospheric exit results in the minimal recircularizing impulse, 


AVj, for single impulse maneuvers. Furthermore, the sensitivity of AV 2 to variations in 
the flight path angle is seen to be high. 

This result stresses the importance of achieving the appropriate conditions at atmo- 
spheric exit. However, off-nominal atmospheric conditions can produce significant pertur- 
bations in the actual trajectory of the vehicle. Significant variations in the atmospheric 
density profile have been observed in shuttle flights. Such variations from the standard at- 
mosphere can result in deviations from the nominal trajectory and perturb the atmospheric 
exit conditions. Since stochastic nonlinear optimization techniques are not currently prac- 
tical, an alternative is to define an optimal control problem which can generate optimal 
trajectories from the current actual state to the desired final state at atmospheric exit, and 
thus adjust to off-nominal atmospheric conditions. This strategy requires fast and reliable 
algorithms for solving two-point-boundary-value problems (TPBVP), which requires fur- 
ther investigation beyond the current study. However, the ability to update the trajectory 
would result in small rather than large variations in the exit parameters. 

In Section 2, a reduced order model of the equations of motion is developed. This 
second order model uses the vehicle’s total energy as the independent variable instead of 
time. Reduction of the order has the advantage that it reduces the order of the TPBVP 
to be solved. Furthermore, it recognizes that coplanar orbit transfer is a problem of 
transferring the vehicle from one energy level to another. The choice of the control value 
is intuitively more a matter of how much energy the vehicle must lose rather than what 
time it is. 

In Section 3, an optimal control problem to transfer the vehicle from an initial state 
to a specified final state is formulated. The necessary conditions are obtained in Section 
4, and optimal trajectories are simulated in Section 5. 
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2 . EQUATIONS OF MOTION 


In this section, we will develop a set of equations describing the motion of the aero- 
assisted orbital vehicle center of mass using total energy as the independent variable rather 
than time. These equations of motion will then be used in the formulation of an optimal 
control problem, where the control will be defined as a function of total energy. As a result, 
although the problem is solved in an open-loop mode, the control is defined in terms of 
energy, so that the control value is selected according to the total energy of the vehicle. 

The equations of motion which will be used describe the dynamics of the c.g. of 
the vehicle in flight within the atmosphere with no propulsive forces used. The general 
equations of motion for this case are given by 


^ _ V cos 7 cos 0 

rcoa<f> 


( 1 ) 


4 > = V cosTfsin ty/r 


( 2 ) 


f = V sin 7 


( 3 ) 


V — 0sin7 + w 2 r cos <^(sin 7 cos 4 > — cos 7 sin <j> sin tf>) 

TTX 


( 4 ) 


. L cos n V g cos 7 . , 

7 = — H h 2u cos^cos 

mV r V 

+(u 2 r/V) cos ^(cos 7 cos <j> -I- sin 7 sin ^ cos iff) 


( 5 ) 
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• Lsmu V 

tp = — — cos 7 cos <p tan <p 

mV cos 7 r 

+ 2w (tan 7 cos sin — sin<£) — (u/ 2 r/Vcos7)sin<£cos<£cosV’ (6) 

where 0 is the longitude, <j> is the latitude, r is the distance from the center of the earth 
to the vehicle center of gravity, V is the velocity with respect to the earth, 7 is the local 
flight path angle, rj) is the track angle, is the bank angle, u> is the singular velocity of the 
earth, m is the mass of the vehicle, D is the drag force and L is the lift force acting on the 
vehicle. 

For the case of coplanar orbital transfer considered here, the general equations of 
motion can be simplified significantly. First, we assume that a lateral regulator control 
system is maintaining the vehicle’s lateral variables near zero by accommodating pertur- 
bations due to atmospheric and other effects. This implies the ability to bank the vehicle 
by small amounts to correct for small perturbations in the lateral variables. Taking the 
initial heading as zero, the latitude, <j>, remains constant and can also be taken as zero. 

The remaining equations of motion are (1), (3), (4) and (5). With the latitude/longitude 
directions as defined above, the motion of the vehicle is along zero-latitude and the posi- 
tion of the vehicle is determined by its longitude and altitude. In the current study, the 
longitude of the vehicle as a function of time is not of interest, and can be eliminated from 
the model since it does not impact the motion. 

Rewriting the remaining equations of motion, we obtain 

r = V sin 7 (7) 
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where the angular rotation of the earth has been neglected and a Newtonian gravitational 
model is used, i.e., 


9 = £ ( 10 ) 

in which the gravitational acceleration, g, is expressed in terms of the inverse square law 
with /2 being the product of the gravitational constant and the mass of the earth. 

Furthermore, we assume a parabolic drag polar form between the coefficients of drag 
and lift: 


C D = C Do + KC 3 l 


( 11 ) 


= = i ClpSV* 


m 


( 12 ) 


*=£ • 


( 13 ) 


Z = l C °' SV2 


( 14 ) 


where p is the atmospheric density, and S is the vehicle’s effective surface area for aero- 
dynamic forces, C Do is the zero-lift or minimum coefficient of drag and K is a coefficient 
principally depending on the vehicle configuration shape. 

Energy as Independent Variable. 

In the model described above, the motion variables r,V and 7 are the dependent 
variables; while the time, t, is the independent variable. However, it is possible to use 
a different independent variable, if some advantages accrue from a different choice. Of 
course, the independent variable cannot be selected arbitrarily. Knowledge of the indepen- 
dent variable must uniquely determine the value of each of the dependent variables. For 



example, altitude cannot be used as the independent variable, since in a standard maneu- 
ver, the vehicle crosses the same altitude twice, once when dipping into the atmosphere 
at a high velocity and again on its way up to the new orbit at a lower velocity. Thus, 
knowing the altitude does not uniquely determine the velocity of the vehicle. 

In the following, we will develop a new model of the motion of the vehicle using its 
total energy, potential plus kinetic, as the independent variable. First note that energy 
is allowable as an independent variable because it is monotonic. Since the aeroassisted 
orbital transfer vehicle (AOTV) will use no propulsive forces during its maneuver inside 
the atmosphere, the vehicle’s total energy monotonically decreases with time due to at- 
mospheric drag. Thus, energy is a one-to-one and invertible function of time; so that the 
vehicle has a given energy level only once during the maneuver. Knowing the vehicle’s 
energy uniquely determines the value of the dependent variables. 

The advantages of using energy as the independent variable are two-fold. First, it 
provides a technique of model order reduction analytically without any approximation. In 
optimal control problems, reducing the model order by one reduces the order of the two- 
point-boundary-value problem which must be solved by two, since the co-state equations 
are also reduced by one. 

The second advantage is that coplanar orbit transfer intuitively is a problem of energy 
management; i.e., it is a problem of going from a high energy level to a lower energy level 
by losing some energy to the atmosphere without overheating. Thus, deciding what control 
value to use is intuitively more a question of how much energy the vehicle has rather than 
what time it is. 

Define the vehicle’s total energy, potential plus kinetic, say E as 

E = \mV 2 — mgr = JmK J — ( 15 ) 

2 2 r 

E-Y1-* 

m 2 r 
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Differentiating E with respect to time and combining with the equations of motion, 


we obtain 


E -V D = -V D . (17) 

As expected, the time rate of change of the vehicle’s energy is work done by the 
aerodynamic forces acting on it; i.e., the force multiplied by the velocity in the direction 
of the force. 

When performing an aeroassisted coplanar orbit transfer, usually the initial energy 
is higher than the final energy. Thus, during the typical maneuver, E progresses in the 
negative direction, whereas the standard variational equations for optimal control problems 
assume the independent variable progressing in the positive direction. To use the standard 
equations, we simply make the change of variables 

e = —E/m , e = -E/m = V D/m (18) 

Now, if the vehicle’s energy and its speed are both known, then its altitude is uniquely 
determined from the potential energy. Thus, it can be found that 



Using the chain rule and (19), the equations of motion (7), (8), (9) can be expressed 
in term of the speed and flight path angle viewed as functions of energy. After some 
manipulation, we obtain 


1 (V 2 + 2e) 3 sin7 

V~ 2A p5V 3 Cp 


( 20 ) 


t _ J_ f £l\ , (K 4 -4c 2 ) cos <7 
7 V* \Cd) 2jipSV 4 C D 


( 21 ) 
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where the prime “ ' ” denotes the derivative with respect to e; i.e., 


de 


( 22 ) 


/ <h 


(23) 


It is seen that the second order model in (20), (21) is sufficient to describe the motion 
of the vehicle. Note that this reduction in the order of the model has been obtained without 
any approximation, by simply using energy as the independent variable. 

The atmospheric density, p, is a function of altitude. For the purpose of the optimiza- 
tion study, the usual exponential form will be used to describe the atmospheric density 
profile with altitude. Thus, 


p{h) = p 0 t h/ & , (24) 

h = r -te , (25) 

where te is the average earth radius and $ the scale height of the exponential atmosphere. 

Since the radial distance, r, is no longer a state variable, the altitude must be expressed 
in terms of the model variables. 


h = 


2A 

V» + 2« * 


(26) 


Thus, the two differential equations (20), (21), the quadratic drag polar (12) and the 
atmospheric density (24), (26) form a reduced order model describing the motion of the 
vehicle center of gravity. 
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3. THE OPTIMIZATION PROBLEM 


The maneuver considered is to transfer the AOTV from high earth orbit to low earth 
orbit by grazing through the atmosphere to lose energy. While in high orbit at a radial 
distance of r\, a tangential retro impulse, AVi, is applied, which puts the vehicle into an 
elliptical orbit with perigee at r p . While ideally, r p would be at the top of the atmosphere, 
say r a , realism considerations require a lower altitude which would ensure that the vehicle 
dips into the atmosphere sufficiently to loose the required energy within a reasonable period 
of time. The vehicle then exits the atmosphere at a lower speed of V* and flight path angle 
*U starting an elliptical orbit with apogee, rj. When the vehicle reaches r 2 , a tangential 
circularizing impulse, AVj, brings it to the desired circular orbit. 

It is well-known that in comparison to the all-propulsive orbit transfer, the aeroassisted 
maneuver requires significantly less fuel to achieve the same orbital transfer [4]. Also note 
that although AVi is the larger impulse, the variability of this impulse is quite small; 
i.e., the variation in AVi which achieves an orbit with a perigee anywhere within the 
atmosphere is quite small. On the other hand, the variation in AVj with variations in 
the exit flight path angle and velocity is significant. Furthermore, the minimum AV 2 
impulse occurs when the exit flight path angle, 7 *, is zero, assuming a single impulse from 
atmospheric exit to the low earth orbit [4]. Therefore, it is highly desirable to achieve the 
needed exit conditions to reduce the variability in the amount of fuel necessary to ensure 
the maneuver. 

The atmospheric portion of the orbital transfer is the part of interest in this study. 
An important consideration during the atmospheric maneuver is to maintain the vehicle’s 
skin temperature at acceptable levels. This is directly related to the heating rate produced 
by the atmospheric conditions. On the other hand, off-nominal atmospheric conditions can 
produce significant perturbations. Although nonlinear stochastic optimization is the most 
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direct approach to the treatment of random atmospheric phenomena, it is beyond the 
available resources of the current study. Instead, we will define a deterministic optimal 
control problem which, if necessary, can be solved on-line to adjust to changing atmospheric 
conditions, so that the guidance trajectory may change with varying atmospheres but still 
achieve the desired exit conditions. Alternatively, the trajectories may be computed off- 
line and stored, although stringent storage requirements would be placed on the on-board 
computer. 

To accommodate the various objectives and constraints discussed, an optimal control 
problem may be posed as follows. Since the exit conditions largely determine the fuel 
requirements, the final flight path angle and the speed are considered to be fixed by the 
low earth orbit radius, r j. On the other hand, the initial conditions are determined by 
the high earth orbit radius, ri. The high entry speed, V*, and flight path angle, 
are also considered to be fixed by the particular maneuver. The main objective during 
the atmospheric maneuver, beyond achieving the atmospheric exit conditions, is to avoid 
overheating the vehicle skin. To a lesser extent, it may be of interest to avoid excessive 
shear stress on the skin of the vehicle. To achieve these objectives, the cost function is 
selected as a linear combination of the square of the heating rate and the drag force acting 
on the skin integrated over the entire atmospheric maneuver. 

The control is the coefficient of lift, Cl. The coefficient of drag, Co, is determined 
by Cl through (12). Both positive and negative values are allowed for Cl. Whether a 
negative value of lift is obtained by a negative pitch angle or by a positive pitch angle 
with the vehicle flying upside down would depend on the vehicle and implementation 
considerations. 

Now, the vehicle speed at atmospheric entry determines the initial energy per mass, 
say *i. Alternatively, the high earth orbit energy diminished by the impulsive energy 
of AVi also determines the energy at atmospheric entry. Similarly, let the energy at 
atmospheric exit correspond to ej. Recall that e is the negative of the energy per unit 
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mass, so that «i < cj. Thus, 


V(«,) = Vi l , -,(«,) =l k 


(27) 


V(«s) = V, , -,(e 2 ) = H 


The heating rate is computed using the expression 


(28) 


Q = AfiV a , A = 3.08 X 10~ 4 . (29) 

A number of different expressions are available for the heating rate which in general, would 
depend on further atmospheric variables. Since only the trends are of interest here, (29) 
is considered sufficient. 

We will use a cost function which is a linear combination of the integrated heating 
rate squared and the drag; i.e., 


J = 



dt , 


(30) 


where t\ and are the initial and final times, respectively. 

Since the independent variable is e, J must be expressed in terms of e. From (18), 


— - 1 
dt V D/m 

Using the chain rule and manipulating, 


(31) 


where 



dt = 



t(V n ,C L )<U 


(32) 


iA^ex 

«1 = — * — » c, = c. 


( 33 ) 
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It is interesting that while the final time, t 2 > in the original problem would be free, in 
the energy model the final energy is fixed. Since the initial and final (or desired) orbits are 
known, the amount of energy which must be expended in the atmosphere is also known, 
although the duration of the atmospheric maneuver is not fixed, and is a part of the 
optimization. 

Thus, the problem of obtaining aeroassisted coplanar orbit transfer trajectories can 
be posed as the optimal control problem of minimizing the cost function, J, in (32) while 
satisfying the constraints (20), (21), (27) and (28). 


12 


4. NECESSARY CONDITIONS 


Using standard variational calculus texts [14], [15], it is possible to determine con- 
ditions which are necessary for optimality. Following this approach, the Hamiltonian, M, 
is 

H{V,~i,C Ly p v ,p n ,e) =g{V,7,C L ) + pvY' + p 7 V , (34) 

where p is the co-state vector defined by 

Pv * * -|^( v »7,C £ ,pv,P 7 ,e) • (35) 

Substituting (32), (20) and (21) into (34), the Hamiltonian for the problem is found 
to be 


H(V,7>C L ,pv,Pi,e) = hi(V,^,p v ,p^,e) — 


h2[V,7,pv,Pi) +MV»7,p v ,P'r) 


(36) 


= *\V 3 - Pv a in If 


[V 2 + 2 e} 2 

2fipSV 3 


+ PACOS'/ 


[V A - 4 e 2 } 
2flpSV 4 


(37) 


^2(V,Tf,pv,p 7 ) = y 


(38) 


hsiV^pvjP-t) = — y PV 


(39) 


The differential equations for the co-state vector can be obtained from (35)-(39). 
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Pv 


o V2 

-3?i — — I- Pv 

C D 


+p 7 


{-^+ 2 (* fl/2 


16e 3 l 

V s 


+ jLViv- 

+ V 3 /*J\2 

cos 7 rv 2 — 2e 


2 upSCo 


+ 


1/2 _ 3e A ain ^ 

V*/*j2iipSC D 

2 cos 7 ^ ?2 


[V 2 + 2eJ j)pSV z C D j V 2 


2 sin 7 \ 

&>SV 2 C D J 

(40) 


. [V 2 + 2e] 2 cos7 [V 4 + 4e 2 ]sin7 

2 jipSV 3 C D 2 ppSV*C D 


(41) 


Note that the atmospheric density, p, is a function of V and e; so that the rate of 
change of p with V directly enters the co-state equation for pv- 


-{ 


;} = 


ifiV 


(42) 


dV x p(V,e)> p(V3 + 2 e) a p(V, e) 

From Pontryagin’s minimum principle, the minimal cost occurs when the Hamiltonian 
is minimized within the allowable control set while the state and co-state are on the 
optimal trajectory; i.e., minimize H{V*,7*,CL,py,p*,e) over the allowable set of CV s. 
First consider the optimal Cl with no constraints on Cl- 


dX 

dC L 


(v~,y ,Ci, p?,,p;,«) = -2 - 2K ^ 


(43) 


where h\ and h 2 are (37) and (38) respectively, evaluated on the optimal trajectory. 
Setting (43) to zero results in 


Cr= jti ± . tax 

c ^ VUa/ 


Cm 

K 


(44) 


Note that the negative sign in front of the discriminant always produces a negative Cl, 
while the positive sign results in a positive Cl value. Observation of (36) shows that the 
minimal X occurs when 


h* 2 C L < 0 


(45) 
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since Cd is always positive. It follows that the root which corresponds to the minimum is 
given by 


• < 46 > 

Now, suppose that Cl is limited to be within [Cimim C^, max ], where Ci, maz is positive 
and CLmin is negative. The minimal Cl for this case is easily obtained by analyzing the 
gradient dM jd Cl in (43). Rewriting this gradient in the form 


d* ( Kh\)C\ + (2 Kh\)C L - h\C Do 

dC L ~ cl * (47) 

Suppose hi is negative, then the gradient has two zeroes, Cl- and Cl+ corresponding 

to the sign selected in (44). Note that Cl- is negative while C&+ is positive. The basic 

shapes of the gradient and the Hamiltonian lead to a value of C£ > 0 which is limited by 

C Lmax- A similar analysis for the case of h\ being positive results in 


Cl> 

CLmin — C[, < C^ maz 

CLmaxi 

Cl > Climax 

CLmin > 

Cl < 


where Cl is given by (46). 


(48) 


The necessary conditions for the optimal trajectory are given by the state equations 
(20), (21), the co-state equations (40), (41), and the control equations (46), (48), with the 
initial and final state satisfying (27) and (28). Thus, the necessary conditions specify the 


two-point-boundary-value problem given above. The sufficiency of these conditions is not 


treated here. 
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5. NUMERICAL RESULTS 


To obtain the optimal guidance trajectories resulting from the optimal control problem 
posed, the necessary conditions were solved using a standard two-point-boundary-value 
problem (TPBVP) solver. The TPBVP solutions were obtained with shooting techniques 
[16] using the OPTSOL software package. 

The problem considered was a typical coplanar orbit transfer from high earth orbit to 
low earth orbit using aeroassist to achieve the maneuver. The initial circular orbit is at an 
altitude of 22,366 km over the earth surface, while the low earth orbit is at an altitude of 
715.6 km. The atmospheric entry conditions resulting from this initial orbit were selected 
to be a speed of 10 km/sec and flight path angle of -6° by choosing a target perigee at a 
radial distance of 6406.5 km from the earth center. 

The atmospheric exit conditions are specified by the speed of 8 km/sec and flight 
path angle of 0.01 rad or 0.57° to ensure a slightly positive flight path angle to exit the 
atmosphere. The top of the atmosphere was selected at 127 km or a radial distance of 
6498 km. The atmospheric scale height was set at 7.5 km with the zero altitude density 
p 0 corresponding to 7.769 x 10 10 kg/km 3 . 

The vehicle parameters Cdo and K were set respective values of 0.05 and 1.4. The 
mavimum lift-to-drag ratio for the vehicle was 1.9. The effective mass to vehicle area ratio 
used was 300 kg/m 3 . 

Since the initial condition of the state, y, is specified, solving the TPBVP consists 
of finding the initial co-state values which will drive the 4 th order state/co-state system 
of differential equations to the desired final state. Since the co-state equations integrated 
forward, are usually unstable, the solution of TPBVP’s is a difficult problem. Significant 
convergence problems were, in fact, encountered in solving the necessary conditions. 
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The sensitivity of the optimal initial co-state vector to changes in the problem pa- 
rameters was found to be high. A consequence of this sensitivity is that the radius of 
convergence of the shooting algorithm for the problem considered was relatively low. So 
that when parameters such as the drag coefficient Cqo or the scale height, of the at- 
mospheric density are varied by small amounts, the algorithm does not converge; this was 
found to be the case particularly at lower values of the scale height. Although a complete 
study was not made, multiple shooting did not significantly modify this situation. On the 
other hand, in many cases, the rate of convergence of the single shooting algorithm was 
fast, requiring under ten iterations. 

A parametric study of the optimal guidance trajectories for different linear combina- 
tions of the heating rate versus the drag terms in the cost function was performed. The 
optimal trajectories obtained are shown in Figures 2-5. The proportions of the heating 
rate (squared) and drag were varied by fixing at .001 while takes on the values of 
2.0, 1.0, 0.6 and 0.0. Recall that when ?2 vanishes the cost function minimizes only the 
heating rate term. As Zi increases, the cost function contains greater proportions of the 
drag force term so that the shear stress on the skin is also included as as objective. 

The basic features of all the trajectories are similar. The speed is reduced from 10 
km/sec to 8 km/sec slowly at first, then at a higher rate until reaching approximately 8.2 
km/sec. At that point, the curve flattens considerably, slowly moving towards its final 
value at 8. Similarly, the flight path angle is increased until it reaches nearly + 1.4°. At 
this point, a rather sharp or decisive reversal of the trend brings the flight path angle to 
a flat curve until atmospheric exit. 

The heating rate increases as the vehicle dips into the atmosphere. However, it remains 
under 100 W/cm 3 throughout the atmospheric maneuver. This level is satisfactory, as 
much higher rates can be accommodated. The coefficient of lift starts near a value of two 
and remains there initially until it drops and settles near a negative value of -0.5. The 
trajectories shown here correspond to the unconstrained control case. Due to convergence 
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problems and time constraints, the constrained control cases were not obtained. 

As the proportion of the heating rate versus drag is increased, the essential character 
of the trajectory remains unchanged. The main difference is seen in the duration of the 
maneuver which increases as the drag term is phased out. Also note that the final flat 
portion of the trajectory is lengthened while the prior portion is slightly shortened in time. 
However, the heating rate is largely unchanged and remains safely under 100 W/cm 2 in 
all the trajectories. The altitude profile tends to become slightly more flat at the end of 
the maneuver when ffj = 0. 
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6. SUMMARY 


The atmospheric portion of the trajectories for the aeroassisted coplanar orbit transfer 
have been investigated. The equations of motion for the problem are expressed using a 
new reduced order model using total vehicle energy, kinetic plus potential, as the indepen- 
dent variable rather than time. The order reduction is achieved analytically without an 
approximation of the vehicle dynamics. 

In this model, the problem of coplanar orbit transfer is seen as one in which a given 
amount of energy must be transferred from the vehicle to the atmosphere during the 
trajectory without overheating the vehicle. An optimal control problem is posed where 
a linear combination of the integrated square of the heating rate and the vehicle drag is 
the cost function to be minimized. The necessary conditions for optimality are obtained: 
These result in a 4 th order two-point-boundary-value problem. 

A parametric study of the optimal guidance trajectory in which the proportion of the 
heating rate term versus the drag varies is made. The problem considers transferring the 
vehicle from an orbit at an altitude of 22,366 km to one at an altitude of 715.6 km in a 
two impulse aeroassisted maneuver. Sensitivity and convergence problems of the shooting 
algorithm are discussed. Simulations of the guidance trajectories are presented. 
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Figure 1. Typical Aeroassisted Coplanar Orbital Transfer 
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Figure 2.a. Optimal Guidance Trajectory for = .001, c 2 = 2.0 
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Figure 2.b. Optimal Guidance Trajectory for Cj = .001, c 2 
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Figure 2.c. Optimal Guidance Trajectory for ?i = .001, c 2 = 2.0 
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Figure 2.d. Optimal Guidance Trajectory for ii = .001, cj = 2.0 
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Figure 2.f. Optimal Guidance Trajectory for c x - .001, c 2 - 2.0 



Figure 2.g. Optimal Guidance Trajectory for c x = .001, c 2 = 2.0 
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Figure 2.h. Optimal Guidance Trajectory for Si — .001, cj — 2.0 
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Figure 2.i. Optimal Guidance Trajectory for £i = .001, = 2.0 
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Figure 2.j. Optimal Guidance Trajectory for i\ = .001, c 2 = 2.0 
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Figure 3.a. Optimal Guidance Trajectory for = .001, e 2 = 1.0 
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Figure 3.b. Optimal Guidance Trajectory for €i = .001, C 2 — 1.0 
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Figure 3.d. Optimal Guidance Trajectory for i\ = .001, 
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Figure 3.f. Optimal Guidance Trajectory for — .001, Cj — 1.0 
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Figure 3.g. Optimal Guidance Trajectory for = .001, e 2 = 1.0 
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Figure 3.h. Optimal Guidance Trajectory for Cj = .001, 



Figure 3.i. Optimal Guidance Trajectory for = .001, e 2 = 10 
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Figure 3 .j. Optimal Guidance Trajectory for Si — .001, £3 = 1.0 
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Figure 4.a. Optimal Guidance Trajectory for i x = .001, e 2 
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Figure 4.b. Optimal Guidance Trajectory for 2i = .001,^2 — 0.6 


130.00 



(^) aaruinv 


45 


Figure 4.c. Optimal Guidance Trajectory for = .001, C 2 = 0.6 
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Figure 4.d. Optimal Guidance Trajectory for l\ — .001,22 = 0.6 
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Figure 
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Figure 4.f. Optimal Guidance Trajectory for Ci = .001, C 2 = 0.6 



Figure 4.g. Optimal Guidance Trajectory for c x - .001, c 2 = 0.6 
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Figure 4.h. Optimal Guidance Trajectory for cj = .001, cj 



Figure 4.i. Optimal Guidance Trajectory for C\ = .001, c 2 = 0.6 
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Figure 4.j. Optimal Guidance Trajectory for C\ = .001, = 0.6 



Figure 5.a. Optimal Guidance Trajectory for c x — .001, c 2 = 0.0 
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Figure 5.b. Optimal Guidance Trajectory for ci = .001, C 2 = 0.0 
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Figure 5.c. Optimal Guidance Trajectory for Ci = .001, C 2 = 0.0 
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Figure 5.d. Optimal Guidance Trajectory for Cj = .001, — 0.0 


200000.00 



AilSNaa *FUV 


57 


Figure 5.e. 




Figure 5.f. Optimal Guidance Trajectory for = .001,52 = 0.0 




Figure 5.g. Optimal Guidance Trajectory for Ci = .001, c 2 = 0.0 
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Figure 5.h. Optimal Guidance Trajectory for c\ = .001, — 0.0 



Figure 5.i. Optimal Guidance Trajectory for Si = .001, cj = 0.0 
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Figure 5.j. Optimal Guidance Trajectory for = .001, c-i = 0.0 
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